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ABSTRACT 

We study the possibility of using the entire probability distribution function (PDF) 
of the aperture mass M ap and its related cumulative probability distribution function 
(CPDF) to obtain meaningful constraints on cosmological parameters. Deriving com¬ 
pletely analytic expressions for the associated covariance matrices, we construct the 
Fisher matrix and use it to estimate the accuracy with which various cosmological 
parameters can be recovered from future surveys using such statistics. This formalism 
also includes the effect of various noises such as intrinsic ellipticity distribution of 
galaxies and finite survey size. The estimation errors are then compared with the ones 
derived from low order moments of the PDF (variance and skewness) to check how 
efficiently the high M ap tail can be used to constrain cosmological parameters such as 
O m , erg and dark energy equation of state w& e . We find that for future surveys such as 
JDEM the full PDF does not bring significant tightening of constraints on cosmology 
beyond what is already achievable by the joint use of second and third order moments. 
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1 INTRODUCTION 

Weak lensing surveys can be used as a very efficient probe 
to constrain the background cosmology as well as contents 
of the Universe (see Munshi & Valageas 2005 and references 
therein). Typically low order quantities such as two- and 
three-point functions or their reduced forms, variance and 
skewness, or the entire PDF are used to obtain cosmolog¬ 
ical constraints from observational data (Bernardeau et al. 
2004). In weak lensing studies, observables such as re, 7 or 
more commonly M ap and their low order moments are ex¬ 
tensively studied (e.g. Jarvis et al. (2004) for recent results 
for M a p). Use of M ap has the additional advantage of be¬ 
ing able to separate gravity induced shear signals (“Elec¬ 
tric” modes) from the (“Magnetic”) modes due to various 
systematics such as point spread function. The one-point 
probability distribution function of M ap encodes informa¬ 
tion regarding non-Gaussianities at all orders (Bernardeau 
& Valageas 2002) and thus can be useful to pinpoint back¬ 
ground cosmological parameters by breaking degeneracies 
which appear at the level of two-point correlation functions 
(Schneider et al. 2002). We develop an analytical formalism 
to study the covariance of binned PDF and employ Fisher 
formalism techniques to study covariance of estimation er¬ 


ror of cosmological parameters from realistic weak lensing 
surveys. The purpose of this study is twofold. Firstly we 
check the estimation error associated with cosmological pa¬ 
rameters while using binned PDF as the primary observable. 
Secondly we use these results to infer how much tightening 
of constraints, if any, can be achieved in general by using 
higher order informations regarding non-Gaussianities. This 
has been a topic of discussion in the recent years using the 
lowest order non-Gaussian statistics, e.g. bi-spectrum or re¬ 
lated collapsed one-point skewness (Takada & Jain 2002, 
Kilbinger & Schneider 2004, Munshi & Valageas 2005). Thus 
our results based on PDF extend such studies to include all¬ 
orders of non-Gaussianity. This letter is organised as follows: 
in section 2 we discuss the analytical results concerning co- 
variance structure of the PDF and CPDF. Borrowing results 
from Valageas et al. (2004) we show how the Fisher matrix 
can be constructed from PDF and CPDF data. In section 3 
the numerical results are provided. Section 4 is left for dis¬ 
cussion of our results. 
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2 ANALYTICAL RESULTS 

2.1 Covariance of PDF of aperture-mass M ap 

We recall in this section the analytical results presented in 
Valageas et al. (2005). The aperture-mass M ap can be writ¬ 
ten as a function of the tangential shear as (Kaiser et al. 
1994; Schneider 1996): 

M ap = J dt? Q Ma p(i?) 7t(t?), ( 1 ) 

with (using the same filter as Schneider 1996): 

.(£)’(.- J). 

Then, in order to measure the aperture-mass A/ ap within a 
single circular held of angular radius 9 „, in which N galaxies 
are observed at positions Dj with tangential cllipticity et.j, 
we can use the estimator Al defined by: 
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Pj ~ ~N~N~ ^ 1 Pj = AjPj = -jy- ^ ' 1 j(n), ( 7 ) 
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where 1 j ( n ) is the characteristic function of the interval Ij = 
[Alj,Mj+i[ of width Aj, applied to the value M(n) of M 
measured in the cell n: 

1 j(n) = 1 if Mj < M(n ) < Mj+ 1 , lj(n) = 0 otherwise.(8) 

Here we restrict ourselves to non-overlapping intervals Ij. 
The estimators Fj are slightly more convenient than the Pj 
as they do not explicitly involve the widths Aj which can 
vary with j or even be infinite (on each side of a binning of 
[A/ m in, A4 max ]). Then, from the sets {Pj} or {Pj} we obtain 
an histogram which provides an approximation to V(AI): 

_ __ r M i +1 

(Fj) = Fj with Pj = / d MV(M). (9) 

J Mj 

Thus, for small enough Aj we have Pj = Fj/Aj ~ P[(A/j + 
Mj +i)/2]. Next, assuming that different cells are well sep¬ 
arated so as to be uncorrelated the covariance matrix C F j F 
of estimators {Pj} is simply: 


In the case of weak lensing the observed complex cllipticity 
e is related to the shear 7 by: e = 7 + e„, where e* is the 
intrinsic cllipticity of the galaxy. Assuming that the intrin¬ 
sic ellipticities of different galaxies are uncorrelated random 
Gaussian variables, the cumulant of order p of Al is: 
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(4) 


where 8 P} 2 is the Kronecker symbol and a 2 = (e 2 ) is the dis¬ 
persion of the intrinsic cllipticity of galaxies. Since the in¬ 
trinsic ellipticities are Gaussian and we neglected any cross¬ 
correlation with the density held they only contribute to 
the variance of the estimator Al (note that M 2 is a biased 
estimator of (M 2 p ) because of this additional term). The 
quantity p measures the relative importance of the galaxy 
intrinsic ellipticities in the signal. They can be neglected if 
p > 1. Any Gaussian white noise associated with the de¬ 
tector can also be incorporated into the expression 0 by 
adding a relevant correction to <r 2 . Finally, from eq. 0 we 
obtain for the generating function </?m of the estimator M: 
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where we defined as usual the characteristic function 1 pm 
of any random variable A I by the logarithm of the Laplace 
transform of its PDF V(M): 


V(M) = 
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Thus, the PDF of the estimator V(AI) is related to the PDF 
of the aperture-mass P(M ap ) through eq.0. Of course, for 
small p we recover the Gaussian (i.e. = —y 2 /2) as Al 

is dominated by the galaxy intrinsic ellipticities, whereas for 
large p we recover y?M ap {y) as A I is dominated by the weak 
lensing signal (AL ~ M ap ). 

Thus, each circular field of angular radius 9 S yields a 
particular value for the quantity M defined in eq.0. If the 
survey contains N c such cells on the sky, we can estimate 
the PDF T(A1) through the estimators Pj or Pj defined by: 


C[j F = (FFj) 

-{ Fi )(Fj)= SijFi N f iFj . 

( 10 ) 

In addition to the Pj which measure the PDF we have also 
considered the estimators Uj (upper) and Lj (lower) associ¬ 
ated with the CPDF: 
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Their covariance matrices are simply: 


/~iTJU _ ^max(ij) UiUj ^LL ^min^j) LiLj ^ 
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The second (product) term is of rank one but the presence of 
the first term makes the matrices non-degenerate. It is also 
possible to construct joint data vectors but this will not be 
pursued here as they all carry the same information as the 
set {Pj}. Note that if the N-p intervals run over ] — oo,oo[ 
so that Pj = 1 the covariance matrices are degenerate 
as there are only N-p — 1 independent variables. Thus, in 
practice we consider the set {Pj} with j = 1,.., Np — 1. In 
this case we also have Lj + Uj = 1 which yields Gjj = C¥j U . 


2.2 Parameter estimation using the noisy PDF 


To obtain the estimation error for cosmological parame¬ 
ters 0 a measured from the observables Pj, Uj or Lj , which 
we generically denote by Xj, we employ the Fisher matrix 
formalism commonly used in the litterature (e.g. Tegmark 
et al. 1997). The information regarding estimation errors 
(AO 2 ) 1 / 2 and their cross-correlations (AQ a AQp) are en¬ 
coded in the inverse of the Fisher matrix T a p which can be 
constructed from the covariance matrix C of the data vector 
and its derivative w.r.t. cosmological parameters 0 a : 
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Figure 1. The covariance structure ( 'f 1 of the estimators {/■)} 
of the non-linear PDF V(M) (left panel) and its Gaussian coun¬ 
terpart (right panel). Non-linearities generate large correlations 
among various bins but the covariance remains diagonally domi¬ 
nated. 



Figure 2. The covariance structure of the estimators {Lj } of the 
non-linear CPDF V{< M) (left panel) and its Gaussian counter¬ 
part (right panel). 


Here we ignored the next order correction terms which in¬ 
volve the derivatives of the covariance matrices (Munshi & 
Valageas 2005) as they are much smaller. The cosmological 
parameters that we consider for extraction from the data 
are as and Wde ■ In this letter we concentrate on one 
angular scale d 3 = 2'. Use of other scales would not change 
much the constraints as nearby scales are highly correlated 
and non-gaussianities measured at very small or very large 
angular scales are too noisy (e.g. Munshi & Valageas 2005). 


3 NUMERICAL RESULTS 

3.1 Covariance matrix and the derivative vector 

For all numerical computations we have adopted a SNAP 
like survey strategy and a concordance ACDM cosmology 
as the fiducial model as in Munshi & Valageas (2005). 

In Fig. 0 we plot the covariance matrix C F j F for {F,} 
using N-p = 91 bins from M = —0.05 up to M = 0.08. 



Figure 3. The PDF V(M) (solid-line) and its Gaussian analog 
(dashed-line) are shown as a function of Mj in the upper-left 
panel along with la error bars. The corresponding derivatives of 
the associated estimators { Fj } of these PDF are also displayed 
as a function of Mj in other panels. We show the derivatives 
w.r.t. Wd e (upper-right panel), to (lower-left panel) and to ag 
(lower-right panel). 


For comparison we plot in the right panel the covariance 
which would be obtained for a Gaussian PDF V(M) with 
the same variance. For both cases there is a very pronounced 
diagonal dominance which makes them well-conditioned. Of 
course, the symmetry with respect to M = 0 is broken when 
gravity-induced non-linearites are included. High M tails 
are generated which follow the positive tails of "P(M ap ) and 
T(M) (Munshi et al. 2004, 2005) whereas a sharper decline 
appears at negative M, see also upper left panel of Fig. [I] 
We plot in Fig. |21 the covariance matrix Cf L for the CPDF 
(here C LL = C uu ). Note that being integrated quantities 
the covariance matrix of {Uj} or {Lj} is more narrowly fo¬ 
cussed along the diagonal but we found that in practice the 
covariance matrix C FF is no more difficult to invert than 
C LL . 

The PDF V(M) is shown in the upper-left panel of 
Fig. EH along with the derivatives of the data vector {Fj } 
with respect to the three cosmological parameters Wde, 
and ag. The absolute value of the derivatives shows a well 
pronounced maxima near M ~ 0 and two wings on each side. 
The symmetry is broken by inclusion of non-Gaussian effects 
as shown by the solid curves. Derivative vectors asymptoti¬ 
cally reach zero for high and low values of M, in bins where 
the PDF becomes small and noise dominated and no useful 
cosmological information can be extracted. We can already 
see that the dependence on Wde is much smaller than the 
sensitivity to S2 m or ag which implies that errorbars for Wde 
will be larger, as will be checked in Fig. 2] below. 


3.2 Parameter constraints using the PDF 

From the the covariance and the derivative of the data vec¬ 
tor {P, } we compute the Fisher matrix to study parame¬ 
ter degeneracies. We consider two parameter pairs, {D m , as} 
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Figure 4. The contribution to the error-ellipses from various seg¬ 
ments of the estimated PDF is shown for the cosmological param¬ 
eter pairs (left panel) and {Q ln , ay } (right panel). The 

dotted line corresponds to the M < 0 sector of the PDF and the 
dashed line to the M > 0 sector. The solid line corresponds to 
the entire range of {Mj}. All other parameters are assumed to be 
perfectly known. All contours in these panels as well as the next 
plots represent 3<r contours. 
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Figure 5. The contribution to the error-ellipses from low-order 
moments are presented for {Dm. (left panel) and (D Trl . a 3 { 

(right panel). The dotted lines represent constraints from variance 
only whereas the dashed lines represent constraints from skewness 
only. The solid lines show the joint constraints. 


where parameters are negatively correlated and {fl m ,uide} 
where they are positively correlated. We also use the two 
lowest order moments (variance and skewness) as data vec¬ 
tors to construct independent Fisher matrices. The latter 
are obtained from the unbiased estimators M v defined as: 

M p = 'y ' Qitap )"-Qm s . p ( 3j p ) etji...ey p (15) 

P Ul,--dp) 

where the sum runs over all sets of p different galaxies 
among the N galaxies enclosed in the angular radius 6 a 
and ( N) p = N(N — 1)...(N — p + 1). This ensures that 
( M p ) = (Ma P ) (contrary to ( M v )). Thus the variance and 
skewness of M ap are not estimated from the PDF V(M). 
The procedure is outlined in Munshi & Valageas (2005) and 
will not be repeated here. As shown by the comparison of 
Figs. m and|S]we find that the constraints from the whole 
PDF are only slightly tighter than the ones constructed from 
these two lowest order moments. The similarity in shape is 
expected since our model for the non-linear PDF is actually 
fully parameterized by the variance and skewness of the un¬ 
derlying density field (see Munshi et al. 2004). Hence the 


n ffi -cr. 
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Figure 6. The estimation error AO Q for two cosmological pa¬ 
rameter pairs are shown, {D m , u; ( j e } in right panel and {f 2 m , a 3 \ 
in left panel (hence this corresponds to la errorbars). The lines 
which reach constant asymptotes for high positive values of Mj 
are for the estimator Lj and the ones which reach constant asymp¬ 
totes for low negative values of Mj correspond to the estimator 
Uj. These asymptotes match exactly with the results based on 
the full estimator set { Fj }, as CPDF and PDF carry equivalent 
information. 


PDF and these two lowest order moments carry the same 
information within our model. The good agreement of our 
analytical predictions with numerical simulations (Munshi 
et al. 2004) shows that this must be true to a large extent 
for any realistic model. In a similar fashion, for the Gaussian 
approximation we noticed that the parameter constraints 
from only variance are in good agreement with the ones de¬ 
rived from PDF (which is then fully defined by its variance). 
However, it was not obvious that the size of the contour area 
would be so close since the low-order estimators {M 2 , M 3 } 
of eq. EH might have been more or less noisy than the data 
vector { F-j }. Indeed, the set {Fj } contains many more points 
than the two-point data vector {M2,M3} but it is based 
on the estimator M which is slightly more affected by the 
galaxy intrinsic cllipticities (e.g. although the variance of 
M 2 is equal to the variance of the estimator M2 the vari¬ 
ance of M 3 is somewhat larger than for M3). In particular, 
we checked that using as data vector {(M 2 ), ( M 3 )} (i.e. the 
two lowest order moments of the PDF V(M)) yields confi¬ 
dence areas which are somewhat larger than for {M2, M3}. 
Thus, it appears that although the full data vector { F 3 } 
is able to yield as good constraints as those obtained from 
{M2, M3} the tails of the PDF are too noisy to provide sig¬ 
nificantly tighter constraints. Note that a drawback of V(M ) 
is that it is more difficult to evaluate in surveys which con¬ 
tain many holes while low order moments may be recovered 
by integrating the measured correlation functions. 

We have also investigated the amount of information 
which can be extracted from the restriction of the PDF to 
M < 0 (dotted lines in Fig .0 or M > 0 (dashed lines). We 
can note that the negative M sector of the PDF generates 
ellipses which are more closely aligned with those from the 
variance analysis alone whereas the positive M sector of the 
PDF produces ellipses which are more closely aligned with 
those from the skewness analysis. This is due to the fact 
that the skewness is mostly sensitive to the large positive 
tail of 'P(Map) which is broader than the negative tail (e.g. 
Munshi et al. 2004). On the other hand, the positive sec¬ 
tor M > 0 alone already provides reasonable constraints on 
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Figure 7. The estimation error A(-),, for two cosmological pa¬ 
rameter pairs are shown for the estimators Fj, {Q rn . w^ e } in right 
panel and {S4 m , erg} in the left panel. A group of six neighbouring 
bins is used for each of these points to reduce the scatter in A(-),,. 


cosmology. We also noticed that the cross-correlation terms 
in the covariance Cf F do not play any significant role and 
for a SNAP like survey can safely be ignored as far as es¬ 
timation of cosmological parameters is concerned. Since the 
CPDF data Uj and Lj are derived from the original data vec¬ 
tor { Fj } they cannot provide additional cosmological con¬ 
straints. However, they can be used as a useful consistency 
check. We have checked that they give exactly the same re¬ 
sults as the PDF data vector { Fj }. 

Finally, we display the estimation error AO a for the 
pairs {fl m , os} and {D m , Wde} associated with different parts 
of the CPDF (Fig. or of the PDF (Fig.Q. The “U-shape” 
in Fig. □is consistent with the upper-left panel of Fig.0 the 
tails of the PDF are too noisy to bring significant cosmolog¬ 
ical information. On the other hand, the small decrease of 
AO q near Mj ~ 0.02 for V(< M) (solid lines in FiglHl is 
related to the secondary peak near Mj ~ 0.02 in the deriva¬ 
tives shown in Fig. □ This feature contains significant cos¬ 
mological information which moreover goes beyond the mere 
variance as its shape differs from the one obtained within a 
Gaussian approximation. 


ing several angular scales. However, as seen in Munshi & 
Valageas (2005) for a space based survey such as SNAP 
with reasonably good sky coverage and high number den¬ 
sity of source galaxies most of the useful information can in 
effect be obtained by studying one particular angular scale 
(e.g. 9 S = 2' in case of SNAP, compare Fig. □with Figs.14, 
15 in Munshi & Valageas 2005). Indeed, nearby scales are 
highly correlated and do not provide additional information 
whereas very large and very small angular scales are more 
affected by noise (due to the intrinsic ellipticity distribution 
of galaxies or the finite size of the survey). 

The results provided here are for a monolithic survey 
strategy but real surveys will have more complicated topol¬ 
ogy. However our results can provide clues as to what extent 
the surveys can probe non-Gaussianity to extract meaning¬ 
ful constraints on cosmological parameters by using not only 
the first few lower order moments but the entire PDF. 

We have not considered a Wiener filter based approach 
to reconstruct the PDF from noisy data as suggested by 
Zhang & Pen (2005) for convergence k maps but such an ap¬ 
proach can very easily be implemented by using the covari¬ 
ance matrices we have constructed. A complete Wiener filter 
based approach using compensated filter A/ ap as presented 
here will be investigated elsewhere. In a different context 
shape of the non-Gaussian PDF was used by Huffenberger & 
Seljak (2005) to separate the kinetic-SZ contributions from 
primordial CMB. Similar approach can be useful in separat¬ 
ing gravity generated E mode using M ap PDF from various 
systematics.Techniques presented here can also be applied to 
the case of weak lensing of diffuse background such as CMB 
(Kesden et al. 2002) and high-redshift 21 centimeter radi¬ 
ation from neutral hydrogen during the era of reionization 
(e.g. Cooray 2004). Results of such analysis will be reported 
elsewhere. 
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4 DISCUSSION 

We have checked the possibility of using the entire Af ap PDF 
to constrain cosmological parameters as opposed to the use 
of its lower order moments which is more prevalent in the lit¬ 
erature. We find that for a SNAP like survey the PDF does 
not yield significantly tighter constrainsts than those derived 
from the variance and skewness alone, despite the low-noise 
space based survey strategy. This also implies that higher 
order moments would not contribute either to further tight¬ 
ening of error ellipses. However, much larger surveys with 
high level of source galaxy distributions might still benefit 
from measurement of kurtosis, but higher number density 
requires increasing the depth of the survey which in turn 
makes the PDF itself more Gaussian. Our results only take 
into account the volume averages of higher order correlation 
functions and did not propose to quantify non-Gaussianity 
beyond collapsed one point objects. It is possible to extend 
our study by taking into account redshift binning or tomog¬ 
raphy which we will report elsewhere, as well as combin¬ 
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